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We show that the universal a-attractor models of inflation can be realized by including an auxiliary 
vector field A^i for the Starobinsky model with the Lagrangian f{R) = R-\- 7?^/(6M^). If the same 
procedure is applied to general modified f{R) theories in which the Ricci scalar R is replaced by 
R-\-A)_iA'^-\-l3XfiA^ with constant /3, we obtain the Brans-Dicke theory with a scalar potential and the 
Brans-Dicke parameter ojbd = /3^/4. We also place observational constraints on inflationary models 
based on auxiliary vector modified f{R) theories from the latest Planck measurements of the Cosmic 
Microwave Background (CMB) anisotropies in both temperature and polarization. In the modified 
Starobinsky model, we find that the parameter /3 is constrained to be /3 < 25 (68 % confidence level) 
from the bounds of the scalar spectral index and the tensor-to-scalar ratio. 


I. INTRODUCTION 

The inflationary paradigm [1, 2] has been the backbone 
of high energy cosmology over the past three decades. For 
the realization of inflation, we require the existence of at 
least one additional degree of freedom to the Einstein- 
Hilbert action. A canonical scalar field with a nearly flat 
potential can play such a role [3, 4]. In modified grav¬ 
itational theories, a scalar degree of freedom generally 
emerges as a result of the breaking of gauge symmetries 
present in General Relativity [5-8]. 

The first model of inflation, which was proposed by 
Starobinsky in 1979 [1], is based on the modification 
of gravity with the Lagrangian f{R) = R-\- R?/{(ilvR), 
where R is the Ricci scalar and M is a constant hav¬ 
ing a dimension of mass. In general, the f{R) gravity is 
equivalent to the Brans-Dicke (BD) theory [9] with the 
BD parameter wbd = 0 [10]. The propagation of a scalar 
degree of freedom in /(i?) gravity is particularly trans¬ 
parent in the Einstein frame where a canonical scalar field 
(j) evolves along a potential of gravitational origin [11]. 

The Starobinsky model gives rise to the Einstein-frame 
potential with a nearly flat region responsible for infla¬ 
tion [6, 12]. In this model the scalar spectral index rig and 
the tensor-to-scalar ratio r of primordial density pertur¬ 
bations generated during inflation are given, respectively, 
by rig — 1 — 2/N and r ~ 12/N‘^, where N is the number 
of e-foldings on scales relevant to the CMB temperature 
anisotropies [13-15]. The Starobinsky model is consis¬ 
tent with the recent joint analysis of the Planck temper¬ 
ature data [16] and the B-mode polarization data from 
BICEP2/Keck array with the Planck maps at higher fre¬ 
quencies [17]. The tensor-to-scalar ratio is constrained 
to be r < 0.08 at 95% confidence level (CL) from such 
a joint analysis. The large-held models like chaotic in- 
hation and natural inflation are now in tension with the 
CMB data [18]. 

Recently, there have been numerous attempts to em¬ 


bed the Starobinsky model in the framework of super¬ 
gravity [19]-[29] or ghost-free higher-derivative gravita¬ 
tional theories [30, 31]. The bottom line is how to 
build up the Einstein-frame potential similar to the form 

= Uo(l - where Mpi = 2.435 x 

10^® GeV is the reduced Planck mass. In a particular 
version of supergravity where the inhaton is a part of a 
vector multiplet it is possible to construct a generalized 
potential of the form U(((>) = Vb(l - e"V2/(3«) 
where the parameter a is inversely proportional to the 
curvature of the inhaton Kahler manifold [22]. This was 
dubbed the a-attractor model in which the inhationary 
period is followed by the reheating stage with the inhaton 
oscillations around ^ = 0. 

In the limit a —>■ c» the potential of the a-attractor 
model is approximately given by V{(j)) oc (p, so it is 
equivalent to that of the quadratic potential in chaotic 
inhation [4]. For 1 < a < oo the tensor-to-scalar ratio is 
in the range 0(10“®) < r < 0(10“®), with rig inside the 
95 % CL observational contour constrained by the Planck 
data [18, 32, 33]. In Ref. [34] the authors derived the 
same potential as that of the a-attractor model by gen¬ 
eralizing the Starobinsky model in the framework of the 
BD theory and they placed observational constraints on 
the model from the WMAP 7yr data. 

In this paper we show that the a-attractor model arises 
by introducing an auxiliary vector field and replac¬ 
ing the Ricci scalar R with R -|- -\- l3'S/fj,A^ in 

the Starobinsky Lagrangian, where is the covariant 
derivative. This procedure can be extended to the gen¬ 
eral f{R) Lagrangian. The resulting theory is equivalent 
to the BD theory characterized by the BD parameter 
wbd = /3^/4 with one scalar propagating degree of free¬ 
dom. 

In light of the recent release of the Planck tempera¬ 
ture and polarization data, we also put observational con¬ 
straints on inflationary models in the framework of auxil- 
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iary vector modified f{R) theories. Our analysis not only 
encompasses the a-attractor model but also the models 
derived by promoting the Lagrangian f{R) = R + cR^ 
(n > 1) to include the auxiliary vector field. This can ac¬ 
commodate a wider class of inflationary models including 
chaotic inflation with the potential oc 

This paper is organized as follows. In Sec. II we re¬ 
view the Starobinsky model and its dual description in 
terms of a scalar degree of freedom ij). In Sec. Ill we 
show how the a-attractor model emerges by modifying 
the Starobinsky model with inclusion of the auxiliary 
vector field. In Sec. IV we extend this prescription to 
general f{R) theories and provide the formulas of infla¬ 
tionary observables associated with the primordial scalar 
and tensor perturbations. In Sec. V we place observa¬ 
tional bounds on the auxiliary vector modified inflation¬ 
ary models from the latest Planck data combined with 
other B-mode polarization data. Sec. VI is devoted to 
conclusions. 


II. STAROBINSKY MODEL AND ITS DUAL 
DESCRIPTION 

The Starobinsky model [I] is described by the action 

pp r 

S=^ J d^Xy/^C{R), ( 2 . 1 ) 

where g is the determinant of the space-time metric 
and C{R) is a function of R of the form 

The discussion given below is already well known in 
the literature [6], but this is useful for the comparison 
with the auxiliary modified Starobinsky model given in 
Sec. III. The model (2.2) possesses an additional scalar 
degree of freedom to that in General Relativity. In or¬ 
der to make this manifest, we consider the following La¬ 
grangian 

It is easy to see that, upon integrating out the field tp, we 
get back to the original Starobinsky model (2.1). Varying 
Eq. (2.3) with respect to F, it follows that 

F = 3M‘^{ip- 1). (2.4) 

Then, the Lagrangian (2.3) can be rewritten as 

C = ipR-^M^{^-l)\ (2.5) 

This is equivalent to the BD theory [9] with the BD 
parameter wbd = 0 and the scalar potential V{(f) = 
{?,/2)M^{(f-l)\ 


In Eq. (2.5) the scalar degree of freedom (f is di¬ 
rectly coupled to the Ricci scalar R. One can trans¬ 
form the action (2.1) with the Lagrangian (2.5) to the 
so-called Einstein frame under the conformal transfor¬ 
mation 'g^i, = Q^{if)gfj_^ [11]. Denoting the quantities in 
the transformed frame as a tilde, we have the following 
relations [6] 

7 =^ = ( 2 . 6 ) 

R = n'^{R + 6au}-6g^^''V^u}V^uj), (2.7) 

where w = InD. We obtain the Einstein-frame action for 

the choice 

n‘^ = (p, ( 2 . 8 ) 

under which the Ricci scalar R does not have a direct 
coupling with (p. Dropping the total derivative term Dw 
in Eq. (2.7) and introducing a scalar field 

0 = yi'ZV/pi Ini^, (2.9) 

the action in the Einstein frame reads 

Ml ~ 1 

( 2 . 10 ) 

where 

V{(l)) = . (2.11) 

Hence the scalar degree of freedom cj), which is the 
gravitational origin, propagates with the kinetic energy 
— (1/2)3^“^The potential (2.11) is sufficiently 
flat for (p larger than the order of Mpi, in which regime 
inflation occurs due to the slow-roll evolution of (j). 

III. AUXILIARY MODIFIED STAROBINSKY 
MODEL 

The auxiliary vector modified Starobinsky model is 
inspired by the supersymmetric extension of Starobin¬ 
sky model [19]-[29]. In the old minimal formulation of 
Af = 1 off-shell supergravity, the Weyl multiplet consists 
of the vielbein e“, the gravitino an auxiliary vector 
Afj_ and an auxiliary complex scalar S. The embedding 
of Starobinsky model in the old minimal supergravity 
is obtained by coupling a chiral multiplet to the Weyl 
multiplet. The supersymmetric Starobinsky model can 
be recast into the form of a scalar-tensor theory by inte¬ 
grating out the auxiliary fields. In particular, integrating 
out the auxiliary vector field generates the kinetic term 
for the imaginary part of the complex scalar in the chiral 
multiplet^. 


^ A detailed performance of this procedure can be found in the 
Sec. 5 of [29]. 
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We apply the similar mechanism here by coupling an 
auxiliary vector field to the Starobinsky model in a spe¬ 
cific way such that the auxiliary vector field does not 
generate new degrees of freedom. As we will show below, 
integrating out the auxiliary vector field modifies the ki¬ 
netic term of inflaton. The resulting theory written in 
the Einstein frame coincides with the a-attractor model 
proposed in Ref. [22]. In analogous to the supersymmet¬ 
ric extension of Starobinsky model, the auxiliary vector 
modified Starobinsky model thus provides a gravitational 
origin for the designed scalar potential in the a-attractor 
model. We can also apply the same mechanism to aux¬ 
iliary vector coupled f{R) theories and obtain a class of 
generalized a-attractor models (see Sec. IV A). 

We would like to stress that our model is inspired by 
the Af = I off-shell supergravity, but it does not di¬ 
rectly come from a supersymmetric scenario with a SUSY 
breaking mechanism. Hence we do not take into account 
the effect of gravitinos for the cosmological dynamics. 
Construction of a SUSY breaking a-attractor model with 
the effect of gravitinos taken into account is beyond the 
scope of our paper. 

We start with the Lagrangian of the form C = R + 
where ,3 is a constant. Note that the 
coefficient in front of the term A^A^ has been fixed to 
1. When higher derivative terms are included, can 
pick up kinetic terms such as and , so 

the auxiliary vector starts to propagate. If we would like 
to allow kinetic terms for A^ but still wish to keep A^ 
as an auxiliary vector in the higher derivative extended 
model, then the action of the auxiliary vector modified 
Starobinsky model has to take the following form: 

f 

s=^ J d‘^xy/^C, (3.1) 

where 

C = i? + A^A^ + 

In order to see that this model gives rise to only one 
scalar degree of freedom, we perform the similar analysis 
to that performed in the previous section. We first write 
Eq. (3.2) as 

F - R-A^Af^ - PV^A^^) . (3.3) 

Varying the Lagrangian (3.3) with respect to A^ and F, 
respectively, it follows that 



(3.4) 

3M^{(p - 1). 

(3.5) 


The equation of motion (3.4) of the auxiliary vector field 
demonstrates that on-shell, the vector field is equivalent 
to the gradient of the scalar field. Therefore, there are 


no dynamical spin-1 degrees of freedom in our model. 
Consequently, the apparent presence of a vector field in 
our model does not spoil the homogeneity and isotropy 
of the universe. 

Substituting the relations (3.4)-(3.5) into Eq. (3.3) and 
dropping a total derivative term, we obtain the (Jordan- 
frame) Lagrangian 

C = (3.6) 

This theory is equivalent to the BD theory with the BD 
parameter 

wbd = , (3.7) 

so the auxiliary vector model (3.2) possesses one scalar 
degree of freedom. 

Under the conformal transformation 
with = ip, the action in the Einstein frame reads 

(3.8) 

where ^ is a canonical scalar field defined by 

^ Yfpi In p . (3.9) 

The potential V{4>) is given by 

V{(j)) = (^1 - , (3.10) 

where 

32 2 

a = 1 J- — = 1 + —wbd • (3A1) 

This is equivalent to the a-attractor model studied in 
Ref. [22]. Setting /3 = 0, we recover the Starobinsky 
model described by the potential (2.11) in the Einstein 
frame. 

Around ^ = 0, the potential (3.10) is approximately 
given by V{<j)) ~ M‘^(fP/(2a). The graceful exit from 
inflation to reheati ng natu rally occurs after the field (j) 
enters the regime ^J2/{?^a)(j)/Mp\ <C 1. In Fig. 1 we plot 
the evolution of the field </> and the Hubble parameter H 
versus the cosmic time t in the Einstein frame for /? = 10. 
During inflation, the field evolves slowly along the poten¬ 
tial (3.10) with a nearly constant Hubble parameter. For 
the parameter /3 and initial conditions chosen in Fig. 1, 
the end of inflation is characterized by the field value 
(j)f ~ O.SMpi with the number of e-foldings N ~ 55. This 
shows good agreement with the analytic estimation given 
in Sec. VA [see Eq. (5.5)]. 

As we see in Fig. 1, the Universe exits from the infla¬ 
tionary epoch to the reheating stage driven by the os¬ 
cillation of (p. The field p exhibits a damped oscillation 
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data combined with other data. 


IV. AUXILIARY VECTOR MODIFIED f{R) 
THEORIES AND INFLATIONARY 
OBSERVABLES 


The discussion in Sec. Ill can be extended to more gen¬ 
eral auxiliary vector modified f{R) theories. In this sec¬ 
tion we shall perform such an analysis and then provide 
the formulas of the primordial power spectra of scalar 
and tensor perturbations generated during inflation. 


A. Auxiliary vector modified f{R) theories 

The auxiliary vector modification to f{R) theories re¬ 
quires replacing R with R + + j3V^A^^. Thus the 

model is given by the action (3.1) with 

C = f{R + A^A^^ + ^V^A^). (4.1) 


Figure 1: Evolution of the scalar field (j) and the Hubble pa¬ 
rameter H in the Einstein frame for /3 = 10. We choose the 
initial conditions 4> = lOMpi and 0 = 0 at t = 0. In this case 
the inflationary period lasts with the number of e-foldings 
about 55, which is followed by the reheating stage with the 
oscillating infiaton field. 


around the potential minimum at ()> = 0. Since the aux¬ 
iliary vector only gives rise to the change of the kinetic 
term of (p in Eq. (3.6), it does not modify the cosmologi¬ 
cal dynamics after the field (j) stabilizes at the minimum 
of the Einstein-frame potential (3.10), i.e., after reheat¬ 
ing. In Appendix A we estimate the time at the onset 
of radiation-dominated era. After this epoch, the energy 
density of radiation dominates over that of </>. 

Since the field 4> behaves as a massive oscillating scalar 
around the potential minimum, the basic mechanism of 
reheating is similar to that in the Starobinsky model 
[6, 35, 36] apart from the fact that the energy scale of 
the potential gets lowered by the factor 1/a. The modi¬ 
fied kinetic term in the Jordan frame can be interpreted 
as the modified shape of the potential in the Einstein 
frame. As we will see in Sec. V A, this modification of the 
Einstein-frame potential gives rise to the change of CMB 
observables relative to those in the Starobinsky model. 
Especially, the larger value of the tensor-to-scalar ratio 
r caused by the modification of the kinetic term of ip is 
a distinguished observational feature to discriminate be¬ 
tween the a-attractor model [0(10“^) < r < 0(10“^)] 
and the Starobinsky model [r = 0(10“^)]. 

The observational constraints on the potential (3.10) 
were discussed in Ref. [34] with the WMAP 7yr data and 
in Ref. [22] with the Planck lyr data. In Sec. V we shall 
place observational bounds on the same model as well as 
more general models from the latest Planck temperature 


We rewrite this Lagrangian of the following form 

C = f{F)-p{F-R-A^A^-py^A^^). (4.2) 

Varying Eq. (4.2) with respect to A^ and E, we obtain 


A, = ^PV,p, 

(4.3) 

/.f(A) = p. 

(4.4) 


Here and in the following, a comma in the lower index de¬ 
notes the partial derivatives with respect to scalar quan¬ 
tities represented in the index, e.g., = df /OF. The 

quantity F depends on p through the relation (4.4). 

Substituting Eq. (4.3) into Eq. (4.2), it follows that 

C = pR- - [pF{p) - f{F{p))] . (4.5) 

This is equivalent to the BD theory with the same BD 
parameter as Eq. (3.7). Compared to Eq. (3.6), the scalar 
potential in the Jordan frame is generalized to Vj(p) = 
pF{p) - f{F{p)). 

Introducing the scalar field (f) as Eq. (3.9) and carrying 
out the conformal transformation 'g^^, = p , we obtain 
the Einstein-frame action (3.8) with the scalar potential 


M^i _ /3r 

V(</) = -^e V 3oTtj;r 


2 <?i> 

F — e 




(4.6) 


<!> 

where f,F{F) = p = “n. The a-attractor model, 

which corresponds to the potential (3.10), is a special 
case of a larger class of auxiliary vector modified f{R) 
theories. When /3 = 0 we have = 0 from Eq. (4.3), so 
that the Lagrangian (4.1) recovers that of f{R) theories. 
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B. Inflationary observables 


Let us study inflation for the theories described by 
the action (3.1) with the Lagrangian (4.1). We assume 
that the background is described by the flat Friedmann- 
Lemaitre-Robertson-Walker metric with the line element 
ds^ = —dt^ + a^{t)5ijdx'^dx^, where the scale factor a{t) 
depends on the cosmic time t. We consider scalar and 
tensor metric perturbations on this background. 

In the Einstein frame the action is given by Eq. (3.8) 
with the potential (4.6). In Refs. [37] it was shown that 
the inflationary observables associated with linear scalar 
and tensor perturbations are invariant under the confor¬ 
mal transformation. 

The spectral index of scalar perturbations with the 
primordial power spectrum Vs is defined byns = I-|- 
dlnVs/dlnk, where fc is a comoving wavenumber. We 
also introduce the tensor-to-scalar ratio, as r = VhlVs, 
where Vh is the primordial power spectrum of tensor per¬ 
turbations. Under the slow-roll approximation during in¬ 
flation, these observables are given by [38] 


where 



* - 127r2M6 U| ’ 

(4.7) 

Us = 1 — 6e + 2rj, 

(4.8) 

r = 16e, 

(4.9) 

V j ’ V ■ 

(4.10) 


As long as the slow-roll condition is satisfied, the analytic 
estimations (4.7)-(4.9) are accurate enough to confront 
inflationary models with the CMB observations. Defin¬ 
ing the tensor spectral index as Uj = dhiVu/dhik, the 
following consistency relation holds [38] 

r = —8nt. (4.11) 


We define the number of e-foldings N = ln[a(t/)/a(t)] 
in the Jordan frame, where a{t) and a(t/) are the scale 
factors at the moments t and tf respectively. The lower 
index “/” represents the values at the end of inflation. 
We identify the held value </)/ by the condition e{(j)f) = 1. 

The number of e-foldings is a frame-independent quan¬ 
tity by properly choosing the observer’s reference frame 
[39]. On using the relations d = Da and dt = iddt 
[11, 12], the Hubble parameters H = {dd/dt)/d and H = 
(da/dt)/a in the two frames are related with each other 
as H = [H + {dn/dt)/id]/n. Since we are considering 
the choice (2.8), the number of e-foldings N = Hdt 
can be expressed as N = Hdi -|- ln((/j/(/j/)^/^. On 
using the slow-roll approximations ~ V and 

3Hd(f>/di ~ in the Einstein frame, it follows that 


The number of e-foldings associated with the CMB tem¬ 
perature anisotropies corresponds to 50 < A < 60 [38]. 
On using Eq. (4.12), the inflationary observables (4.8) 
and (4.9) can be known as functions of N. 


V. OBSERVATIONAL CONSTRAINTS FROM 
THE LATEST CMB DATA 


We put observational constraints on several different 
inflationary models that belong to the class of auxiliary 
vector modified f{R) theories. 

We employ the bounds in the {ns,r) plane derived by 
the latest Planck CMB temperature data (Temperature- 
Temperature (TT), Temperature-E-mode (TE), E-mode- 
E-mode (EE) correlation power spectra) and a first re¬ 
lease of the B-mode polarization data [16]. The Planck 
mission also performed the joint analysis by taking into 
account the B-mode maps from BICEP2 and Keck Ar¬ 
ray with the Planck maps (BKP) at higher frequencies 
where the emission of dust dominates [17]. This study 
showed that there is no statistical significant evidence 
for the detection of primordial gravitational waves. Still, 
the latest BKP analysis put tighter upper bounds on the 
tensor-to-scalar ratio r than those derived by the Planck 
data alone [18]. 

The likelihood analysis of the Planck mission is 
based on expansions of the scalar and power spectra 
of the forms Vs{k) = 

and Vh(k) = At(/c//c*)"’*+(“*/^)respectively, 
where Og,* = dns^t/dlnk are the runnings of the spec¬ 
tral indices and fc* is the pivot wavenumber. Since there 
is no significant evidence for the large deviation of as^t 
from 0, the standard slow-roll prediction of inflationary 
observables is consistent with the CMB data. 

In Fig. 2 we show the 68 % CL and 95 % CL observa¬ 
tional contours in the (ns,r) plane derived by the latest 
Planck temperature data as well as the BKP and Baryon 
Acoustic Oscillations (BAO) data. From the Planck TT, 
TE, EE and low-multipole temperature polarization data 
(denoted as “lowP”), the tensor-to-scalar ratio r is con¬ 
strained to be r < 0.15 (95% CL) [18]. Combination 
of the BKP cross-correlation with the Planck TT-t-lowP 
data gives a tighter bound r < 0.08 (95 %CL). Inclusion 
of the BAO data leads to the shift of Us toward larger 
values (as in the figure 1 of Ref. [33]). In the follow¬ 
ing we shall place observational constraints on concrete 
auxiliary vector modified /(i?) models. 


A. Auxiliary vector modified Starobinsky model 


Let us begin with the model (3.2), i.e., 

fiF) - F + g^2 ’ 


(5.1) 


N = 


V 






d(j) + 


1 

^/^ Mpi 


(4.12) 


in the Lagrangian (4.2). Since the potential in the Ein¬ 
stein frame is given by Eq. (3.10), the observables (4.7)- 
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Figure 2: The 68 % CL (inside) and 95 % CL (outside) obser¬ 
vational contours in the {ns,r) plane derived from the joint 
data analysis of Planck TT-|-lowP-|-BKP-|-BAO (thick solid) 
and Planck TT, TE, EE-|-lowP (thick dashed). The pivot 
scale is chosen to be fc, = 0.002 Mpc“^. We also show the 
theoretical curves of the auxiliary vector modified Starobin- 
sky model (3.2) as functions of /3 for N = 50 and 60 (thin 
solid). 


(4.9) reduce to 


mass scale M is constrained to be 

for 0 = 0(1), (5.7) 

In the presence of the coupling /3, both M and r are larger 
than those in the the Starobinsky f{R) model (a = 1). 

In the limit that cr ^ 1, inflation occurs in the re¬ 
gion around a; = 1, so the potential (3.10) reduces to 
V{4)) ~ (fp'/{2a). This means that, for increasing 
/3, the observables (5.3) and (5.4) approach the val¬ 
ues of the quadratic potential, i.e., Us — 1 — 2/N and 
r ~ S/N. Since the scalar power spectrum is given by 
Vs — A^^M^/(67r^Q(Mpj), the Planck normalization gives 

~ 6.6 X for a » 1. (5.8) 

Mpi V 55 / 

In this regime the mass scale M is higher than that for 
a = 0(1). 

In Fig. 2 we plot the theoretical curves in the {ns,r) 
plane as functions of /? (ranging 0 < /3 < 10®) for N = 50 
and 60. The quadratic potential is outside the 95 %CL 
observational contours. For = 60 the joint data analy¬ 
sis of Planck TT-t-lowP-fBKP-l-BAO gives the following 
bounds 


/3< 25 (68% CL), (5.9) 

/3< 66 (95% CL). (5.10) 

For N = 50 the Starobinsky model (/3 = 0) is out¬ 
side the 68 % CL region (mainly due to inclusion of the 
BAO data), but it is still inside the 68% CL contour con¬ 
strained by Planck TT, TE, EE-l-lowP. 


3a / A/ y (l-x)^ 
1287r2 \MpJ ’ 

(5.2) 

^ 8a;(a;-|-l) 

3a(l — a;)2 

(5.3) 

64a;^ 

(5.4) 

3a(l — x)2 ’ 


where x = e ^ number of e-foldings 

(4.12) reads 



where Xf = {2y/^ — 3q:)/(4 — 3q:). 

When the parameter a = 1 + /3^/6 is of the order of 1, 
the inflationary epoch corresponds to the regime in which 
(j) is larger than Mpi, i.e., a; <C 1. Since N ~ ia/{4,x) in 
this case, it follows that 


* 247r2aMpj ’ 




(5.6) 


From the Planck normalization Vs — 2.2 x 10 ® [16], the 


B. Power-law model 

We proceed to the power-law model given by 

/(F) = , (5.11) 

where m and n are positive constants. Since F"“^ = 
ip/{nmP'P~'^^) from Eq. (4.4), the Einstein-frame poten¬ 
tial (4.6) reduces to 

V{(t>) = 

. . 

The positivity of the potential requires the condition n > 
1. The power-law inflation (a cx V with p > 1) [40] can 


be realized for 3a > [{n — 2)/{n — 1)]^. 
have 

In this case we 

n-1 

" 3a(n-l)2’ 

(5.13) 

r = 8(1 - Us). 

(5.14) 


The Harrison-Zeldovich spectrum corresponds to the 
limit n —>■ 2 or cx —> oo. 


n — 1 


2n”/(' 


— Mpim exp 


n — 2 


























7 



s 


Figure 3: The theoretical curves of the model (5.15) with 
= 60 in the (us,r) plane for /I = 0,10,10“^. These curves 
are plotted as as functions of the power n (> 1). The theo¬ 
retical prediction for the n = 2 case is shown as a thin dotted 
curve. For smaller n, the tensor-to-scalar ratio gets larger. 
The 68 % CL and 95 % CL observational contours are the same 
as those plotted in Fig. 2. 


The theoretical values of Ug and r are on the line (5.14) 
in the (jis^r) plane. This line is outside the 95 %CL 
regions shown in Fig. 2, so the power-law model (5.11) is 
disfavored from the data. 


C. Generalization of the auxiliary vector modified 
Starobinsky model 


Finally we study the following model 

/(F) = (5.15) 


where m and n 2) are positive constants. In this case 
the potential in the Einstein frame is given by 


n 


V{cj^) = - 1 ) 

2n”/("“i) P V / 

(5.16) 

Since we consider inflation in the regime ((> > 0, the pos¬ 
itivity of the potential requires that n > 1. 

For given n and /?, we numerically compute the field 
value (/)/ at the end of inflation according to the condi¬ 
tion e(0/) = 1. From Eq. (4.12) we identify the field 
value corresponding to = 60 and then evaluate the 
observables (4.8) and (4.9). 

Let us first discuss the case /3 = 0. Since inflation 
occurs in the regime where (j) is bigger than Mpi, the 


large deviation of the power n from 2 spoils the flat¬ 
ness of the potential. In Fig. 3 we plot the theoretical 
curve in the (rig, r) plane for iV = 60 as a function of n. 
The Starobinsky model (n = 2) is shown as a black cir¬ 
cle. For smaller n the tensor-to-scalar ratio gets larger, 
whereas the scalar spectral index reaches a maximum 
value rig = 0.991 around n = 1.92 and then turns into 
decrease. When n > 2, the inflationary observables are 
particularly sensitive to the deviation from n = 2 because 
of the appearance of the potential maximum in the Ein¬ 
stein frame. Erom the Planck TT-t-lowP-|-BKP-|-BAO 
joint analysis we obtain the following bound 


1.980 < n < 2.015 (/? = 0), (5.17) 

at 95 % CL. Hence only the tiny deviation from n = 2 
is allowed for the consistency with the CMB data^. The 
mass scale m constrained by the Planck normalization 
is similar to that given in Eq. (5.7), with the correspon¬ 
dence m = y/6M and a = 1. 

In Fig. 3 we also plot the theoretical curve for ,5 = 10 
as a function of n. The qualitative behavior of rig and 
r with respect to the change of n is similar to that for 
/3 = 0. The Planck TT-flowP-|-BKP-|-BAO joint analysis 
gives the bound 

1.75<n<2.39 (/3 = 10), (5.18) 


at 95% CL. The wider range of n is allowed relative to 
the case /3 = 0. This reflects the fact that, for larger 
/3, inflation can occur in the regime where the quantity 
j. — g-\/2/(3“)<^/.^pi is not very much smaller than 1. 
When [3 = 10, for example, the order of x satisfying 
the bound (5.18) is typically 0.1 at = 60, whereas, for 
/3 = 0, X is of the order of 0.01 for n ranging in Eq. (5.17). 

In the limit /3 —>■ oo the epoch of inflation corresponds 
to the regime a: ~ 1 and hence the potential (5.16) can 
be approximated as 


V(cjr) 


n-1 / 2 
2nP \3a 


p/2 


Ml^m^ 



where 


n 



(5.19) 

(5.20) 


This is equivalent to chaotic inflation with the power-law 
potential E(())) oc (jf, so the inflationary observables are 
estimated as 


Vg 

rig 


r 


Cn f m \ 

127r2p2 \MpiJ 


P + 


2N 


4,p 

77’ 


2 



(5.21) 

(5.22) 

(5.23) 


^ This is also related to the fact that the corrections like XnR^ 
(n > 2) to the Starobinsky model f{R) = R-\- need 

to be strongly suppressed during inflation [41]. 


















where c„ = (n—l)[2/(3Q;)]^/^/(2n^’). In the limits n—>• 1 
and n —>■ oo we have p —>■ oo and p —?► 1, respectively. On 
using the relation 4>/Mp\ ~ the Planck normal¬ 

ization constrains the mass scale m, as 

^ ~ 5.1 X (2pAf)-(P+2)/4 . (5.24) 

Mpi 

In Fig. 3 we plot the theoretical curve for (5 = 10'^ 
in the range 1.2 < n < 100. The values of Ug and r are 
very close to those estimated from Eqs. (5.22) and (5.23). 
Provided that n > 3.5, the models are inside the 95 % CL 
region constrained by the Planck TT-t-lowP-fBKP-|-BAO 
data. However, even the linear potential V{(j)) oc (j) (i-6., 
n —>■ oo) is marginally inside the 95 % CL contour, so the 
models with /3 3> 1 are not generally favored from the 
CMB data. 


VI. CONCLUSIONS 

In this paper we showed that the auxiliary vector mod¬ 
ification to the Starobinsky model derived by replacing 
R with a + Afj,A^ + /XVfj,A^ gives rise to the universal a- 
attractor model proposed in the context of supergravity. 
Applying the same prescription to general /(i?) theories, 
the resulting action is equivalent to that of BD theories 
with the BD parameter wbd = /5^/4. Under the con¬ 
formal transformation to the Einstein frame, it is clear 
that one scalar degree of freedom (a canonical field (/>) 
propagates along the scalar potential. 

For the potential with a sufficiently flat region the 
scalar degree of freedom cj) not only leads to inflation at 
the background level, but also the field perturbation Scj) 
can be the source for primordial density perturbations 
relevant to the CMB temperature anisotropies. Using 
the invariance of scalar/tensor perturbations under the 
conformal transformation, the inflationary observables in 
auxiliary vector modified f{R) theories are simply given 
by Eqs. (4.7)-(4.9) with the number of e-foldings (4.12). 

In light of the recent release of the Planck temper¬ 
ature and polarization data, we placed observational 
constraints on inflationary models in the framework of 
auxiliary vector modified f{R) theories. We studied 
three different models: (i) f{F) = F + F^/{6M^), (ii) 
f{F) = m^A-n)pn^ = F + 

(n ^ 2), where F = R + A^A^^ + 

The model (i) is equivalent to the a-attractor model 
with the correspondence a = 1 -I- /3^/6, which recovers 
the Starobinsky model for j3 = Q. From the joint data 
analysis of Planck TT-|-lowP-t-BKP+BAO the parame¬ 
ter P is constrained to be /3 < 25 (68 % CL) for V = 60 
(see Fig. 2). The model (ii) gives rise to the exponential 
potential in the Einstein frame, in which case the the¬ 
oretical line in the {ug , r) plane is outside the 95 % CL 
observational contours. 

For the model (iii) with /3 = 0, the power n is con¬ 
strained to be in the narrow range around n = 2, 
i.e., 1.980 < n < 2.015 (95% CL) from the Planck 


TT-|-lowP-|-BKP-|-BAO data. With increasing /?, the al¬ 
lowed region of n tends to be wider because inflation 
occurs for x = not very much smaller 

than 1. In the limit /3 —>• oo the theoretical values of Ug 
and r are the same as those in chaotic inflation with the 
potential oc which case the model is 

marginally inside the 95 % CL observational contour for 
n > 3.5. 

The issue of super-symmetrization of the auxiliary vec¬ 
tor modified Starobinsky model would be interesting. 
However this is not an easy task, so we leave it for future 
investigation. From the observational side, the possible 
detection of primordial gravitational waves will be able 
to clarify whether or not the Starobinsky model and the 
auxiliary vector modified f{R) models are observation- 
ally favored. We hope that we can approach the origin 
of inflation in the foreseeable future. 
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Appendix A: The onset of radiation-dominated era 

We estimate the time G at which the energy density 
of radiation dominates over that of the field (p for the 
auxiliary modified Starobinsky model. In the Starobin¬ 
sky model {a = 1), this issue was already addressed 
in Refs. [6, 36]. During the oscillating stage of in- 
flaton the potential (3.10) is approximately given by 
U(^) ~ m^0^/2, where m = Hence the dis¬ 

cussion for the model a 1 is analogous to that given in 
Refs. [6, 36] after the replacement of M with m. In what 
follows we estimate the time R briefly. 

To study the particle production during reheating, let 
us consider a massless canonical scalar field y in the Jor¬ 
dan frame. We express the quantum field y in terms of 
the Heisenberg representation: 

(Al) 

where dfc and a\ are annihilation and creation operators, 
respectively. The rescaled field Uk = axk obeys the equa¬ 
tion of motion 

+ k^Uk = U[rpuk , (A2) 

where U{r]) = a^R/6, and rj = J a~^dt is the conformal 
time. The time-dependent term on the r.h.s. of Eq. (A2) 
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leads to the production of x particles with the initial vac¬ 
uum state described by the solution /V^- 

The energy density pr of the held x is associated with 
the Bogoliubov coefficient /?fc = — ^ 

Pr = ^T^k'^dk ■ k\l3k\^, where g* = 

0 (100) is the number of relativistic degrees of free¬ 
dom. During reheating the Ricci scalar evolves as i? = 
0 (1) ^^ sin[m(t — tos)] in the regime m{t — tos) 1, 
where tos is the time at the onset of reheating. Taking 
the time average of oscillations of R, the energy density 
of created particles can be estimated as [6, 36] 


The radiation density (A3) decreases slowly relative to 
(which is proportional to the Held density p^). The 
onset of radiation-dominated epoch (t = tr) is identihed 
by the condition = p^, i.e., 




OS 


4 Mp\ _ 3/2 -^pl 

3C m3 “ 3C“ M3 ■ 


(A5) 


Using the observational constraint (5.8), which is valid 
in the regime a ^ 1, we obtain 


Pr=C-^, (A3) 

where C is a coefficient of the order of 1. The scale factor 
evolves as a oc {t — during the oscillating phase of 

(j), so the evolution of the Hubble parameter squared is 
given by 




4 

9(t - tos)^ ' 


(A4) 


^ 5 X IQi® [nV 1 


(A6) 


Substituting the values C = 0(1) and iV = 55 it follows 
that tr — tos ~ 10“^^ sec. For t > tr the inhaton energy 
density p^ becomes negligible relative to pr, so it does 
not affect the thermal history of the Universe after the 
onset of the radiation era. 
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